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ABSTRACT 

We present high spectral resolution mid-IR observations of SiS towards the C-rich 
AGB star IRC+10216 carried out with the Texas Echelon-cross-Echelle Spectro¬ 
graph mounted on the NASA Infrared Telescope Facility. We have identified 204 
ro-vibrational lines of 28 Si 32 S, 26 of 29 Si 32 S, 20 of 28 Si 34 S, and 15 of 30 Si 32 S in the 
frequency range 720 — 790 cm -1 . These lines belong to bands v = 1 — 0, 2 — 1, 3 — 2, 
4 — 3, and 5 — 4, and involve rotational levels with J ]ow < 90. About 30 per cent of these 
lines are unblended or weakly blended and can be partially or entirely fitted with a 
code developed to model the mid-IR emission of a spherically symmetric circumstellar 
envelope composed of expanding gas and dust. The observed lines trace the envelope 
at distances to the star < 351?* (~ 0" 7). The fits are compatible with an expansion ve¬ 
locity of 1 + 2.5(r/l?* — 1) km s _1 between 1 and 51?*, 11 km s _1 between 5 and 201?*, 
and 14.5 km s -1 outwards. The derived abundance profile of 28 Si 32 S with respect to 
H 2 is 4.9 x 10 -6 between the stellar photosphere and 51?*, decreasing linearly down to 
1.6 x 10“ 6 at 201?* and to 1.3 x 10 -6 at 501?*. 28 Si 32 S seems to be rotationally under 
LTE in the region of the envelope probed with our observations and vibrationally out 
of LTE in most of it. There is a red-shifted emission excess in the 28 Si 32 S lines of band 
v = 1 — 0 that cannot be found in the lines of bands v = 2 — 1, 3 — 2, 4 — 3, and 5 — 4. 
This excess could be explained by an enhancement of the vibrational temperature 
around 201?* behind the star. The derived isotopic ratios 28 Si/ 29 Si, and 32 S/ 34 S are 
17 and 14, compatible with previous estimates. 

Key words: line: identification line: profiles — stars: AGB and post-AGB 
circumstellar matter — stars: individual: IRC+10216 — infrared: stars 


1 INTRODUCTION 

IRC+10216 (CW Leo) is the closes t asymptotic giant 
bran ch (AGB) star to Earth (~ 120 pc; iGroeneweeen et al.l 
2012). During the last decades this star was usually con¬ 
sidered as an i solated star with an effe ctive temperature 
of ~ 2300 K (iRidewav fe Keadvl 1 1988h and surrounded 
by an optically thick C-rich circumstellar envelope (CSE) 
composed of molecular gas and dust, which results from 
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the ejectio n of stellar material at a rate ~ 1 — 5 x 
10~ 5 (e.g., Knapp fe Morrisl 1+85|_ Cernicharo et al.l Il996l : 


ISchoier fc Olofssonl 1200 ll : iDe Beck et al.l 120121 1 . Although 
there were indirect eviden ces suggesting that it could be 
actually a binary system dMauron &; Hugginsl 1999^ 2000|; 
Fong et al.ll2013l: iDecin et al.ll201ll . 120141 : ICernicharo et al l 


2015al c), Kim et alJ 1 20141) recently reported the discovery 


of what could be a faint companion of the AGB star respon¬ 
sible of the ejected matter. 


ICernicharo et al.l 

Although most 
of them form in the outer envelope, the so-called par¬ 
ent molecules (e.g., CO, C 2 H 2 , HCN, SiS, SiO, CS, SiC 2 , 


, Kawaguchi et al. 


200sl Patel et al. 

2011]) 
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Si 2 C) are formed i n warm and dense she l ls near the 
stella r photosphere llKeadv fc Ridgwavl 1 19931 : iBovle et al.l 
119941 : ICernicharo et al.l l201f)l . l2015bl ) where strong dy- 
namical mec hanisms probably exist llHinkle et al] 1 19821 : 


! I Sowe ii 19881 and an active chemistry has been inferred 
ll Willacy fc Cherchnefll 1 19981 : lAgiindez fc Cernicharol l200fJ : 
ICherchnefli 2006l l. The dust formation and the subsequent 


gas a cceleration involved in the development of the envelope 
(e.g.. lGilmanlll972l : lMorrisl[l987l : iBoweijfl 9881 : iFonfrfa et al.1 
I2008TI are supposed to be intimately related to the re¬ 
fractory nature and abundance of the parent molecules, 
in particular to C- and Si-bearing spe cies such as SiS, 
SiC 2 , and Si 2 C (ICernicharo et alji2015bl) . These molecules 
are believed to be highly depleted due to dust through¬ 
out the first 20 — 30i?* outwards from the stellar p hoto¬ 
sphere (|Turnei 19871: Bieeing fc Neuven-Quang-Rieul 


sphere (turner 198/; cicgmg &: INguycn-Uuang-Kicu 
Bieging fc TafallJ 1993l~ Bovle et ahj 19941; iLucas et al.1 


1989 


1995 


Agundez et al.ll2012: iFonfria et alJl20l4l. 

The emission of SiS has been analysed i n the mid- 
infra re d dBovle e t al.| 1993), the s u bmilli metre llPatel et al.l 
rl. 


120091. 120111: iDecin et al.l l20ld. l2014h. a nd the mil- 


limetre frequency ranges (Morris 1975|; 

Henkel et al. 

19831, 19851: Turner 1987; Kahane et al. 

1988, 

200d; 

Bieging & Nguyen-Quang-Rieu 1989; Carlstrom et al. 

199C(; 

Bieging & Tafalla 1993; Lucas et al. 1995); Lucas 

1997j; 

Cernicharo et ah] 2000l, 2014; Fonfria ExDosito et al.1 

200® 

Schoier et al.l 2007!; Agundez et al. 20121: Fonfria et al. 

20 li 

Velilla Prieto et al. 20151). Its abundance throughout the 


middle and outer envelope is reasonably well estimated 
(~ (1 — 2) x 10 -6 outwards from ~ 10—1577* ~ 0.2 —0.3 arc- 
sec) but its value in the innermost envelope is still contro¬ 
versial due to discrepancies of up to one order of magnitude 
between the results of several works. 

The main goal of the current work is to derive a re¬ 
liable SiS abundance profile for the innermost envelope 
of IRC+10216 in order to settle the controversy through 
the analysis of high spectral resolution molecular observa¬ 
tions in the mid-IR, which has been proven as a power- 


single setting covered a spectral range of ~ 0.25 /jm. The 
entire range (11.6 — 13.9 /im) was covered with 10 separate 
settings. We nodded IRC+10216 off the slit for sky subtrac¬ 
tion. The weather was good enough during the observations 
that the median background level measured in the off po¬ 
sitions showed no significant variations. We used a black 
body-sky difference spectrum to correct for the atmosphere. 
The data were re duced with the standard TEXES pipeline 
dLacv et alJl200j ). We normalised each spectrum by remov¬ 
ing the baseline with a polynomial fit. The rms of the ran¬ 
dom noise is ~ 0.15 per cent of the continuum emission. 

Although most of the lines in the spectrum are produced 
by C 2 H 2 , HCN, and their main isotopologues (IFonfrfa et al.l 
l2008h . we have identified 204 lines of 28 Si 32 S, 26 of 29 Si 32 S, 
20 of 28 Si 34 S, and 15 of 30 Si 32 S in the frequency range be¬ 
tween 720 to 790 cm -1 . A hundred and fifteen of the 28 Si 32 S 
lines belong to band v = 1 — 0, 54 to band v = 2 — 1, 13 to 
band v = 3 — 2, 16 of band v = 4 — 3, and 6 of band v = 5 — 4. 
Seventy three of them are unblended or partially blended 
and can be fitted (Fig. [I]). The identified lines of 29 Si 32 S, 
28 Si 34 S, and 30 Si 32 S pertain to band v = 1 — 0. Three, two, 
and one of them, respectively, are unblended or partially 
blended and can be reasonably fitted. The lines of bands 
involving vibrational states with higher energy are below 
the detection limit. No lines of the other SiS isotopologues, 
which are at least 100 times less abundant th an 28 Si 32 S (e.g., 
iKahane et al.l [20001 : ICernicharo et al.l I 2 OO 0 I ) . have been de¬ 
tected in our spectrum. 

We have used the low spectral resoluti on ISO / S WS ob ¬ 
serva tions carried out on 1996 May 31 (ICernicharo et ahl 
Il999ll to determine the properties of the dusty envelope. We 
estimate the uncertainty of these observations due to noise 
and the calibration process in ~ 10 per cent. 


(Keadv et al. 1988 

: iKeadv & Ridgway 19931: Bovle et al.l 

199"3: Fonfria et al.1 

20081). 3 


I 11 this Paper we present the identification of 265 lines 
of the vibrational bands v = 1 — 0, 2 — 1, 3 — 2, 4 — 3, and 
5 — 4 of the silicon sulphide isotopologues 28 Si 32 S, 29 Si 32 S, 
28 Si 34 S, and 30 Si 32 S observed with high-spectral resolution 
(~ 3 — 4 km s _1 ) towards IRC+10216. About 30 per cent of 
them have enough quality to be properly fitted and analysed. 
The observing procedure is included in Section [2] Comments 
on the fitting process of the molecular lines can be found 
in Section [3] The description of the observed SiS lines, the 
results of the fitting procedure, and their analysis are ad¬ 
dressed in Section[]] Finally, a brief summary of the current 
work and our conclusions can be found in Section [5] 


2 OBSERVATIONS 

We observed IRC+10216 with the Texas Echel on-cross- 
Echelle Spectrograph (TEXES; lLacv et al.l l2002l l at the 
NASA Infrared Telescope Facility (IRTF) on 2002 Dec 12. 
We used the TEXES high resolution echelon grating with 
a first order grating as the cross-disperser. In this mode, a 


MODEL OF ENVELOPE AND FITTING 
PROCEDURE 

The observed SiS lines w ere fitted and analyse d with the aid 
of the code developed by IFonfrfa et al.l 120081 ). which solves 
numerically the radiation transfer equation of a spherically 
symmetric envelope composed of gas and dust in radial ex¬ 
pansion. It was successfully used to analyse the mid-IR spec¬ 
tra of C 2 H 2 , HCN, C 4 H 2 , and C 6 H 2 at 8 and 13 /im to¬ 
ward^ IRC+1021+ and the protoplanetary nebula CRL618 
(IFonfrfa et al.ll2008l . l201ll ). The code does not solve the sta¬ 
tistical equilibrium equations, which would require a very 
large amount of CPU time to determine the population of 
all the ro-vibrational levels involved in the calculations. In¬ 
stead, it makes use of rotational and vibrational temperature 
profiles provided by the user. The code solves the radiation 
transfer equation and produces the spectrum of the envelope 
convolved with a given Ga ussian beam. The ne w 3D version 
of this code, presented in IFonfrfa et alJ 1 2(11 ll ). allowed us 
to calculate as well the emission of a simple asymmetrical 
envelope. A deeper insight on how this code works in its 
ID and 3D v e rsions and its benchmarking c an be found in 
IFonfrfa et all J2008I) and IFonfrfa et all J2014I) . 



















































































































































SiS in the innermost envelope of IRC+10216 3 
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Figure 1. 28 Si 32 S, 29 Si 32 S, 28 Si 34 S, and 30 Si 32 S lines in the observed spectrum (black hi stogram). Most of the lines surrounding those of 
the SiS isotopologues belong to C 2 H 2 , HCN, and their isotopologues jFonfrfa et al. 2008). The synthetic spectrum is plotted in red. The 
abscissa axis represents the Doppler velocity ( v = Vlsr — ^sys, Vsys = —26.5 km s ^ ICernicharo et al .|[2000l ; iHe et al.|[2008l ; |Patel et al.l 
12011 I) . The vertical dark grey dotted and dashed lines indicate the velocity of the star (0 km s 1 ) and the maximum and minimum gas 
expansion velocity projected on to the line-of-sight (+U 00 and — Uoo), respectively (Voo = 14.5 km s —1 ). The frequency range shaded 
in light grey in every inset indicates the part of the lines used in the fits. The intensity of some lines with respect to the baseline was 
multiplied by 3 and 5 to improve their visibility (it is indicated at the upper right corner of the corresponding insets as x3 or x5). 
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Table 1. Stellar and envelope parameters involved in the fits 


Parameter 

Units 

Value 

Error 

Ref. 

D 

pc 

123 

— 

2 

a* 

arcsec 

0.02 

— 

3 

-R* 

cm 

3.7 x 10 13 

— 


M 

Mq yr" 1 

2.1 x 10“ 5 

— 

4 

T+ 

K 

2330 

— 

3 

^He 


0.2 

— 

5 

Rin 

R* 

5 

— 

4 

Rout 

R* 

20 

— 

4 

r eX p(l R* ^ r < Rin)* 

km s _1 1 

+ 2.5(r/R* - 1) 

— 

1 

^exp(-^in ^ F < -^out) 

km s —1 

11.0 

— 

4 

^exp^ ^ Rout ) 

km s -1 

14.5 

— 

4 

T k (lR* ^ r < 9R*) 

K 

T*(R*/r)°- 58 

— 

6 

T k (9R* 7 r < 65R*) 

K 

T*(9R*/r) 040 

— 

6 

T k (r ^ 65R*) 

K 

T*(65 R*/r) 12 

— 

6 

Ar(lR* ^ r < Ri n ) 

km s -1 

5 (R*/r) 

— 

7 

A v(r 7 Rin) 

km s _1 

1 

— 

8 

/ac 

per cent 

95 

— 

4 

/SiC 

per cent 

5 

— 

4 

T dust(ll H 


0.70 

+0.11/ - 0.09 

1 

^dust(^in) 

K 

825 

+130/ - 160 

1 

'Tdust 


0.39 

+0.04/ - 0.06 

1 


D: distance to the star; a*: angular stellar radius; R*: linear 
stellar radius; M: mass-loss rate; T*: stellar effective temperature; 
Xft e : abundance of He; i?; n and R out : positions of the inner 
and outer acceleration shells; ?; eX p: gas expansion velocity; T k : 
gas kinetic temperature; Av: line width; /x: fraction of dust 
grains composed of material X, i.e., amorphous carbon (AC) 
or silicon carbon (SiC); r f i usti : dust optical depth along the 
line-of-sight; Tp UHt : temperature of dust grains; 7d u st : exponent 
of the decreasing dust temperature power-law. The parameters 
for which no uncertainty has been provided (—) were assumed 
as fixed throughout the whole fitting process. *The fit to the 
observed SiS lines was better with this expansion velocity field in 
Region I and it was assumed during the whol e fitting process (see 
Sectio n [4.1 .If. References: (1) This wor k (2) iGroenewegen et al] 
l|2012 h f.'j I iRidtrw av fc Ready] 1 1 988|) 14 i [Fonfrfa^^^^^004_(5) 
[Cox ~|200(j) (6) DtM3eck_et al.1 d2012li (7) lAgundez"et~al.l ll20l3 ) 
('8') iHuggins &: Healvl lll98^V 


3.1 Adopted envelope model and spectroscopic 
constants 


In order to analyse the SiS lines identified in the observed 
spectrum, we adopted the model of envelope composed of 
several concent ric shells with a con stant expansion veloc¬ 
ity derived by iFonfrfa et alj (120081) from the analysis of 
the mid-IR spectra of C 2 H 2 and HCN, which are pro¬ 
duced in the same region of the envelope of IRC+10216 
that the observed SiS lines trace. Despite the fact that 
some authors used a different model with a more contin¬ 
uous expansion velocity profile (e.g.. ISc hoier e t al.l [2006b; 
IDecin et al.l I2OI0I . 12014 ; iDe Beck et al. 20121) . the adopted 
model seems to be more suitable to describe the inner layers 
of the envelope llK eadv et al.lll988 l:lKeadv fc Rid awavll 1 99.4; 


Bovle et al .1 1 1994^ Fonfri^etaL I2008L 201^ Cernicharo et all 


201ll . 120131 : lAgundez et al.ll2012l : iDaniel et al.ll2012l ). With 

the aim of supporting our choice, we briefly discuss about 
the quality of the fits using th e sele cted velocity profile and 
that recently derived by IDecin et al] (I2014T) in Appendix 0 
Thus, we assumed that the CSE is divided into three Re¬ 
gions (I, II, and III) separated by two acceleration shells 


located at 5 and 20f?*, which were initially 0.5R* thick each 
(Table [T]). The gas expansion velocity, v e x P , is initially as¬ 
sumed to be 5, 11, and 14.5 km s~ x in Regions I, II, and 
III, respectively. In order to improve the expansion veloc¬ 
ity field close to the star, different thicks of the inner ac¬ 
celeration shell (ending at 5f?*) and expansion velocities 
in Region I were tried. The mass conservation results in a 
gas density profile following very approximately the law oc 
r 2 VgJp, since a negligible amount of the total ejected mass 
(< 0.25 per cent) is lock ed in the dust grains (IDe Beck et al.l 
120121 ; IDecin et al.l l2014 ). We assumed that t here is no atomic 
H in gas phase ( Kna in fe Bowersl 1 19831 : iGlasseoldl 1 19961 : 
lAgundez fc Cernicharo 20061 ). The line width at the stel¬ 
lar surface is assumed to be 5 km s -1 , decreasing to the 
terminal value of 1 km s -1 at 5Jth following the power- 
law oc 1/r (e . g., Huggins^JIeal^ 19861: Keadv et al] 1 19881 : 


ISchoier et al.l 2006al : Agundez et al. 20121) . This line width 
profile is com patible with t he tur bule nt velocity > 3 km s -1 
suggested by Keadv et al.l dl988l) and iMonnier et al.l (120001) 
and the line width of ~ 5 — 8 km s _1 inferred from the vi¬ 
brational excited lines of SiS, CS, and HCN observed in the 
mm and sub-mm wavelength ran ges by IPatel et al.1 (120091 . 
Im3) and lCernicharo et all (l201ll ). which are supposed to be 
formed in the inner layers of the envelope. The kinetic tem¬ 
perature was adopted to follow the prof ile derived fro m the 
analysis of the CO emission bv IDe Beck et ahl (120121) . i.e., 
T k oc r -0,58 between the stellar surface and 9J?*, oc r -0 ’ 4 
from 9.R* to 65J?*, and oc r -1 ’ 2 outwards. The rotational 
and vibrational temperature profiles were assumed to fol¬ 
low r 1 , a dependence on the distance to the star tight ly 
related to the cooling proc ess of the envelope llMcKee et al.l 
119821 ; iDotv fc Leunall997l) . 7 can be different for each re¬ 
gion and for the rotational and vibrational temperatures. 
7 cannot be determined for both the rotational and vibra¬ 
tional temperatures in Region III due to the low sensitivity 
of the model to the emission coming from th at zone. Thus, 
we adopted 7 = 0.55 an d 1.0, respectively llFonfrfa et al.l 
120081 : lAgundez et al .11 20121 ). We assumed SiS is in rotational 
LTE at the stellar photosphere. No dust was considered to 
exist in Region I. Throughout the rest of the envelope, we 
adopted dust grains composed of 95 per cent of amorphous 
carbon (AC) and 5 per c ent SiC, wi th a size of 0.1 /im, and 
a density of 2.5 g cm -3 llFonfrfa et al.ll200sl) . 

The rest frequencies of the observed lines and the en¬ 
ergy of the vibrational states involved in the calculations 
(v ^ 9; Eyjb 5, 9400 K) we re taken from the MADEX 
code dCernicharo et al. 20121) and the fits presented by 
lYelilla Prieto et al.l 1 20151) . The optical constants needed to 
calculate the opac i ty of the dust grains were taken from 
Rouleau & Martin ( 199 ll) for amorphous carbon and from 


Mutschke et al] ( 19991 1 for SiC. 


3.2 Fitting procedure 


The fitting procedure is based on the minimisation of the 
reduced y 2 function defined as 


2 


Xred 


1 


n — p 




, (i) 


where n and p are the number of frequency channels and 
of parameters involved in the fitting process, (F v /F co nt)obs 
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Table 2. Parameters involved in the fits of 28 Si 32 S lines 


Parameter 

Units 

Value 

Error 

Ref. 

ai(lR* < r < R in ) 

xicr 6 

4.9 

+ 0 . 8 / - 0.6 

1 

x(r = Rout) 

xlO -6 

1.6 

+0.3/ - 0.3 

1 

x(r is 50R*) 

xlO -6 

1.3 

— 

2 

Trot (R* ) 

K 

2330 

— 


F r ot (Rin) 

K 

860 

+130/ - 120 

1 

Trot (Rout) 

K 

400 

+70/ - 40 

1 

7rot 


0.55 

— 

2 

T v i bi „ = i_o(R*) 

K 

1560 

+140/ - 120 

1 

T v ib,i>=l-o(Rin) 

K 

613 

+ 11 / - 11 

1 

T v ib, v= 1-0 (Rout ) 

K 

227 

+ 12 / - 12 

1 

Tvi b ,u=2-1 (R* ) 

K 

1290 

+130/ - 80 

1 

T v ib, ij= 2-1 (Rin) 

K 

490 

+18/ - 10 

1 

T v ib, d= 2-1 (Rout ) 

K 

190 

+60/ - 70 

1 

T v i b ,n=3-2 (R* ) 

K 

760 

+80/ - 100 

1 

Tvib,u=3-2 (Rin) 

K 

418 

+ 20 / - 20 

1 

T v ib,t;=3-2 (Rout ) 

K 

153* 

— 

1 

Tvi b ,u=4-3 (R* ) 

K 

550 

+220/ - 230 

1 

T v ib, n =4-3 ( Rin ) 

K 

313* 

— 

1 

T v ib,n=4-3 (Rout) 

K 

117* 

— 

1 

Tvi b ,u=5-4 (R* ) 

K 

312+ 


1 

T v ib,u=5-4 (Rin) 

K 

224* 

— 

1 

T v ib, d= 5-4 (Rout) 

K 

00 

O 

— 

1 

7vi b 


1.0 

— 

3 


x : abundance relative to H 2 ; T ro t: rotational temperature; 7 ro t: 
exponent of the rotational temperature power-law outwards from 
Rout; T v it>: vibrational temperature; 7 v ib : exponent of the vibra¬ 
tional temperature power-law outwards from R 0 ut for all the vi¬ 
brational bands. T ro t and T v ^ are oc r~ 7 , where 7 depends on 
the Region of the CSE and the temperature. The parameters 
for which no uncertainty has been provided (—) were assumed 
as fixed throughout the whole fitting process. * Linearly extrap¬ 
olated from the vibrational temperature of other bands at the 
same position in the envelope. I This parameter is free but the 
fitted lines ar e too weak to de r ive a relia ble error. References : (1) 
This work (2) lAgundez et al.l (120121) (3) iFonfrfa et alj i 20081) . 

and (-F^/Fcont)synth are the flux normalised to the contin¬ 
uum for the synthetic and observed spectra, respectively, 
and (Ti is the uncertainty of the observed normalised flux for 
each channel. In this work, cr* was considered equal to the 
noise rms for all the frequency channels. We assumed a num¬ 
ber of physical parameters related to the envelope model as 
fixed during the fitting process (gas expansion velocity field, 
H 2 density, kinetic temperature, or line width, among oth¬ 
ers; see Tables \l\ and n. The parameters related to 28 Si 32 S 
derived from the fits were the abundance with respect to 
H 2 , the vibrational temperature at the stellar photosphere, 
and the rotational and vibrational temperatures at the in¬ 
ner and outer acceleration shells. The 28 Si/ 29 Si, 28 Si/ 30 Si, 
and 32 S/ 34 S isotopic ratios were estimated from the lines of 
29 Si 32 S, 29 Si 32 S, and 28 Si 34 S as scaling factors of the abun¬ 
dance profile of 28 Si 32 S ('Section 14.211 . The fits to the lines 
of the SiS isotopologues were performed assuming the ro¬ 
tational and vibrational temperatures derived from the fits 
to the lines of 28 Si 32 S. Several other parameters related to 
the dust component of the envelope, such as the dust op¬ 
tical depth or the temperature of the dust grains, were left 
as free parameters in order to fit the continuum and kept 
fixed when the SiS lines were fitted. The fits were carried 
out using all the observed lines from each isotopologue si- 



Figure 2. Dependence of the X 2 e a function around its minimum 
(Xi e d m i n — 60.8) on a given free parameter P. P lmrl is the value of 
the parameter P derived from the best fit. The parameter x(R t < 
r < R n ) is plotted in solid black (•), x(r = -Rout) in dashed black 
(°), Trot (Riii) in solid cyan (*), Trot (Rout) in dashed cyan (*), 
T v ib,u = 1 - o{R*) in solid red (■), T vibi „ = i _ o(Rin) in dashed 
red (□), T vibi „ = i _ 0 (Rout) in dotted red (o), T vibi „ = 2 - l (R*) 
in solid green (a), T v ; b „ = 2 — l (Rin) in dashed green (a), 
T v i b ,„ = 2 - l(Rout) in dotted green (v), T vihiV - 3 _ 2 (R*) in 
solid blue (x), T vib „ _ 3 _ 2 (Rin) in dashed blue (x), and 
T v ; b „ = 4 _ 3 (R*) in dotted blue (+). See the caption of Tables [T] 
and [ 2 ] for the definitions of these parameters. 

multaneously. A number of observed SiS lines were slightly 
offset along the frequency and normalised flux level axes to 
get the best fit. These offsets in the normalised flux level 
axis are included to remove artefacts introduced during the 
baseline removal or due to blends with other lines coming 
from IRC+10216 or weak telluric lines. The offsets in the fre¬ 
quency axis are always smaller than the spectral resolution 
of the observations (~ 3 — 4 km s _1 ) but larger than typical 
with TEXES observations. The offsets have larger ampli¬ 
tude, although still generally less than 2 km s _1 , toward 
higher frequencies. The line list we used is based on exper¬ 
imental measurements and has small errors in our spectral 
range (< 0.5 km s _1 ), as expected for a simple molecule like 
SiS. We see slightly small amplitude frequency offsets when 
we looked at the atmospheric features, suggesting that the 
determination of the frequency scale is less precise at higher 
frequencies for these data. 

The model was sensitive to changes in most of the pa¬ 
rameters related to SiS at distances < 35.R* and to dust 
up to ~ 400.R*. The uncertainties of the free parameters 
(Tables |T] and [ 2 j were calculated by varying them until the 
the difference between the observed and the computed syn¬ 
thetic spectrum was larger than the detection limit of 3<r in 
each channel. This conservative decision directly affects the 
derived uncertainties, which could be overestimated. Dur¬ 
ing this process, all free parameters were varied at the same 
time to account for any possible influence between them that 
could substantially affect the derived uncertainties. These 
errors reflected the sensitivity of the model to each parame¬ 
ter, as can be seen in the shape of the Xred function (Fig. [2]), 
which shows the largest curvature around its minimum (best 
fit) for the most influential parameters on the model (e.g., 
the vibrational temperature of bands v = 1 — 0 and 2 — 1 at 
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the end of the inner acceleration shell, i?.j n ). The minimum 
Xred achieved during the best fit is ~ 60. This large value 
is mostly related to contamination of the SiS lines by other 
features in the frequency intervals used to calculate the Xied 
function. Although these differences do not severely affect 
the global minimisation process due to the large number 
of lines involved, they increase the minimum value of Xred- 
This value could be diminished by modifying the frequency 
interval where the x?ed function is calculated for each line 
to remove what seems a contamination after the end of the 
fitting process, but this choice could compromise the objec¬ 
tivity of the minimisation procedure and was discarded. 


4 RESULTS AND ANALYSIS 
4.1 28 Si 32 S 

Most of the observed lines in bands v = 1 — 0 and 2 — 1 show 
shapes compatible with P-Cygni profiles (Fig. 0. The emis¬ 
sion component cannot be clearly recognised in the profiles 
of the lines of band v = 3 — 2, although this fact can be con¬ 
sequence of strong overlaps with close lines. This component 
is insignificant in the profiles of the lines of bands v = 4 — 3 
and 5 — 4. The lines of band v = 1 — 0 display typical absorp¬ 
tion depths of 10 — 20 per cent of the continuum emission 
while the lines of bands v = 2 — 1, 3 — 2, 4 — 3, and 5 — 4 
are weaker with absorption depths smaller than 5, 2.5, 1.3, 
and 0.7 per cent of the continuum, respectively. The maxi¬ 
mum absorption for all the lines occurs at velocities between 
— 14.3 and —4.5 km s . The absorption component of the 
lines of band v = 1 — 0 with Jj ow < 12 and > 23 show main 
contributions peaking at ~ —13.5 and —9 km s -1 in aver¬ 
age, respectively. The absorption component of these lines 
have widths of ~ 7 km s -1 , even for the lines involving ro- 
vibrational levels with the highest energies. The absorption 
of the lines with 12 < Ji ow < 23 show both absorption con¬ 
tributions (at —13.5 and —9 km s -1 ) and, consequently, the 
largest absorption widths (~ 9 — 10 km s^ 1 ). A width of 
7 km s ^ 1 is also found in the absorption components of the 
lines of band v = 2 — 1 , which display their maximum be¬ 
tween — 10 and —5 km s -1 . The absorption of the lines of 
bands v = 3 — 2, 4 — 3, and 5 — 4 show smaller widths of 
~ 2.5 — 6.5 km s -1 , peaking between —6 and —4 km s -1 . 

The velocity of the maximum absorption and the shape 
of the absorption component of the lines are consequences 
of the expansion velocity of the emitting gas. The width of 
the absorption components can be produced by expansion 
velocity gradients, line broadening due to turbulence, the 
varying thermal line width due to the large temperature gra¬ 
dient, and the different directions along what the absorbing 
gas expands (more important for the gas close to the star). 
Thus, the absorption of the low excitation lines is mostly 
formed far from the star where the gas is colder and it ex¬ 
pands at large velocities. The high excitation lines require 
warm gas and/or an efficient radiative pumping mechanism, 
conditions that are fulfilled close to the star. Hence, the ab¬ 
sorption of the observed high excitation lines indicate that 
the gas in the innermost envelope expands at low velocities. 
The observed medium excitation lines, mostly produced be¬ 
tween the regions where the low and high excitation lines are 
formed, display large absorption widths which suggest the 



Figure 3. Comparison between the observed 28 Si 32 S line 1 — 
0R(17) (black histogram) and the synthetic line calculated with 
three different models depending on the expansion velocity field 
in Region I (blue and red), and on the symmetry of the envelope 
(green). All these lines have been taken from the best global fit 
to the observations. The line in blue has been calculated with a 
constant expansion velocity of 5 km s 1 in Region I. The line 
in red is produced by a spherical symmetric envelope with a gas 
expansion velocity of 1 + 2.5(r/R* — 1) between 1 R* and 5R* (Re¬ 
gion I). The line in green is the synthetic line calculated assuming 
the envelope model of the red line but with an asymmetric vibra¬ 
tional temperature distribution (see Section 14.1.41 . The curves 
in grey depart 3cr from the observations, where a is the rms of 
the noise of the spectrum (— 0.15 per cent of the continuum). 
The widths included in the Figure stand for the FWHM of the 
emission component of the profiles. 


existence of an expansion velocity gradient due to the ab¬ 
sorption of gas at different velocities. The detailed analysis of 
the absorption components of the observed lines performed 
in the current work suggests that the actual expansion veloc¬ 
ity profile and our choice (Section 13.111 are essentially com¬ 
patible but it can be improved close to the star. 


4.1.1 Improving the gas expansion velocity field 

The part of the synthetic SiS line profiles arising from the 
gas in Regions II and III is in a reasonably good agree¬ 
ment with the observations. These calculations assumed the 
expansion velocity field initially considered in our model 
(Section 13.11) . Although a constant gas expansion velocity 
of 5 km s _1 can be adopted as the average velocity in Re¬ 
gion I ( 1 .R* < r < 5-R*), a linear dependence between the 
stellar atmosphere and 5 R* ranging from 1 to 11 km s , i.e., 
Wex P (l R* ^ r < 5.R*) = 1 + 2.5(r/R* — 1) km s _1 , describes 
better the emission component of the 28 Si 32 S lines of bands 
v = 1 — 0 and 2 — 1 (Fig. 0. This velocity gradient gives rise 
to a wider emission component of the low excitation lines, 
although it is worth noting that the observed lines of band 
v = 1 — 0 with J < 40 show a red-shifted emission excess 
with respect to the synthetic lines that is not accounted by 
the improved expansion velocity field (see Figs. 0 and 0 
and Section 0T74] for a discussion about this topic). Lines of 
bands v = 3 — 2, 4 — 3, and 5 — 4 are reasonably well repro¬ 
duced with both expansion velocity fi elds since these lines 
are formed close to the photosphere (lAgundez et al.ll 2012 l : 
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Table 3. Abundance of 28 Si 32 S relative to H 2 


R* 

Ain 

-^out 

> 50.R* Error 

IR Phase 

Ref. 

4.9 

4.9 

1.6 

1.3 

0.08 

1 

3.0 

3.0 

2.4 

1.3 1.5 - 2.0 

Multiepoch 

2 

4.0 

4.0 

4.0 

4.0 4 

0.04 

3 

19 

17 

1.7 

1.7 2-3 

0.05 - 0.24 

4 

51 

33 

2.5 

st 0.9 2 

0.5 

5 

7.5 

7.5 

7.5 

7.5 2-3 

0.09 - 0.29 

6 

The abundances are 

multiplied by 10 -6 . 

Their errors are pro- 

vided 

as multiplicative/dividing factors. 

The pulsation 

phase 

was calculated following the analysis of the infrared light 

curves 

performed bv Monnier et al.l (19981). *See Table [2l (1) This work 

(2) lAgundez 

et al. ( 2012) (3) Decin et al. (20101) (4) Schoier et al. 

(2007 

) f5) Bovle et al. (19941) 16) Bieging & Nguyen-Ouang-Rieu 


19891) 


IVelilla Prieto et al.l 12015). where the expansion velocity is 
expected to be low. 

Hence, we used the improved expansion velocity field 
with the gradient in Region I in our fits. 


4-1-2 8 Si 32 S abundance with respect to H 2 


The best fit to the whole set of lines derived assuming the 
envelope model presented above was achieved with a con¬ 
stant 28 Si 32 S abundance with respect to H 2 by number of 
4.9 x 10 -6 in Region I followed by a linear decrease in Re¬ 
gion II down to 1.6 x 10 -6 at the outer acceleration shell, 
and a further decrease outwards down to 1.3 x 1CP 6 at 
50.R*. The 28 Si 32 S abundance beyond t he outer acce l eratio n 
shell is compatible with the results by Boyle et al.l (1 1994h , 
ISchoier et al.l (l2007l l , and lAgundez_et_ah 


(2012) within a fac¬ 


tor of 2 and with the results by Decin et al. ( 20101 ') within up 
to a factor of 3, a disagreement smaller than the uncertainty 
they estimated for their own results (Table Q3)- However, our 
28 Si 32 S abundance is about 10 and 4 time s sm aller in Region 
I than the results by I Boyle et alj (1 19941 4 and ISchoier et al.l 
(I2007T) . respectively. The large discrepancy of one order of 
magnitude found betwe en our derive d abundance in Region 
I and that proposed bv lBovle et al.l (1 19941 ') is due to signif- 
ic ant differences in t he spectra. In the spectrum presented 
bv lBovle et al.l J 19941 ) only the blue-shifted part of the lines, 
corresponding to the absorption depths, were totally covered 
during the observations while the emission components were 
largely (but not completely) missed. The emission compo¬ 
nent of most of the lines they observed was probably blended 
with telluric lines or features of other molecules existing in 
the envelope of IRC+10216. In addition, it is usually diffi¬ 
cult to determine the baseline in spectral regions with many 
mol ecular lines, such as the region containing SiS (see figs. 1- 
6 in lFonfrfa et ahlfeoPSl ). Thus, the line profiles were assumed 
to display no emission resulting in lines with absorption com¬ 
ponents deeper than they actually were due to the remnants 
of the emission components. The modelling of these lines re¬ 
quired an exceedingly low vibrational temperature and a sig¬ 
nificant increase of the estimated 28 Si 32 S abundance in Re- 
gion I and pa rt of Region II compared to the actual profile. 
ISchoier et alJ (120071 ) derived the abundan ce by fitting previ¬ 
ous ob servations acquired with BIMA by Biegin g fe Tafallal 
(i 19931) . These observations seem not to be affected by any 


instrumental effect suggesting t hat the d isagr eemen t be¬ 


tween our results and those by ISchoier et al.l (1200 7h may 
be produced by an actual variation in the 28 Si 32 S emission. 
Such variation would not be related to the stellar pulsation 
phase, known to influence the emission of sever al molecules 
(ICarlstrom et al] 1 19901 : ICernicharo et all 120141') , since both 
observati ons were carried out close to maximum light (see 
Table [3] Monmer_e£_al. 19981). Th e abundanc es proposed 
by [B ieging fe Nguven-Quang-Ried (1 19891) and iDecin et al.l 


(I2010D . although based on observations acquired near max¬ 
imum light, are affected by large errors that prevent us 
from deriving reliable conclusions about the time evolution. 
Therefore, the observed variation in the 28 Si 32 S abundance 
could be the co nsequence of a larger tim e-scale change in the 
mass-loss rate (ICernicharo et aITl2015a| [cl) or the gas phase 
chemistry close to the star, considering that our observa- 
tio ns were acquired about 1 4 yr after the observations taken 
bv lBieging fe Tafallal (Il993l) . and the ejected gas is expected 
to spend about 10 — 15 yr travelling across Region I. 

The 28 Si 32 S abundance with respect to H 2 cal¬ 
culated with chemical models under thermodynamical 
equilibrium (TE) is larger than ~ 1 0~ 5 in Regi on I 

(I Willacy fe Cherchnefll Il99l : lAgundez fe Cernichard 120061 ). 
i.e., a factor of 2 larger than our results. The TE and chem¬ 
ical kinetic models give constant abundances in most of Re¬ 
gion I and beyond, meaning that gas-phase chemical reac¬ 
tions involving SiS are only expected to p la y an important 
role in the vicinity of the star (e.g., lAgundez fe Cernicharol 
120061 ). Hence, the large decrease in the SiS abundance out¬ 
wards from the inner acceleration zone is probably due 
to depletion on to du st g ra ins, as was suggeste d previ ¬ 
ously ('T urned 1987 : Bieging fe Nguven- Quang-Ri eulll989l; 


Bieging fe Tafalla 1993 : Bovle et al.1 Il994l ; ISchoier et all 


2007 : Agundez et al.1 20121 + 


4-1.3 Rotational and vibrational temperatures 

The rotational temperature profile in the ground vibrational 
state derived from our fits throughout the dust formation 
zone (1R+ < r < 30i?*; e.g.. iFonfrfa et al]l2008l : IDecin et al.l 
I2010I) are similar to the kinetic temperature proposed in 
previous works within the estimated uncertainties (Fig. [I]), 
demonstrating that the 28 Si 32 S molecules in the vibrational 
ground state is under rotational LTE throughout the first 
~ 357?* of the envelope. Moreover, our results suggest that 
the rotational temperature for the vibrational states v = 1 to 
5 can be assumed to be the same for the vibrational ground 
state in the region of the envelope where these vibrational 
states are significantly populated, i.e., at distances to the 
star < 10-R*. 

The vibrational temperatures derived from bands v = 
1 — 0, 2 — 1, 3 — 2, 4 — 3, and 5 — 4 are significantly smaller 
than the kinetic temperature (Fig. [J), which indicates that 
28 Si 32 S is vibrationally out of LTE and its vibrational states 
are, mainly, radiatively populated. Using the gas density 
profile presented in Section l3.il we estimate the gas density 
at the stellar photosphere, n gas {R+), to be ~ 4.2x 10 10 cm -3 . 
The critical density, n cr i t , at this position can be estimated 
with the A-Einstein coefficients of the observed 28 Si 32 S ro- 
vibrational tran sitions, which range from 0.8 to 6.2 s -1 (the 
MADEX code; ICernicharo et ai1l2012l ; IVelilla Prieto et al.l 
2015). and the de-excitation ro-vibrational state-to-state 
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Figure 4. Kinetic, rotational, and vibrational temperatures of 
28 Si 32 S from the stellar photosphere up to 35R*. The kinetic 
temperatures, T^, derived in previous works (references 1, 2, 3, 
and 4) are plotted in solid, short dashed, dashed-dot, and long 
dashed, black lines. The rotational temperature, T ro t, derived in 
the frame of the current work is plotted in solid red. The vi¬ 
brational temperatures, T v ib, derived from ro-vibrational bands 
v = 1 — 0, 2 — 1, 3 — 2, 4 — 3, and 5 — 4 are plotted in blue, ma¬ 
genta, green, cyan, and grey solid lines, respectively. The dashed 
lines in these colours represent possible dependences of the vi¬ 
brational temperatures on the distance to the star b etween the 
stellar photosphere and ~ 1.2R*, where [Cernicharo et alJ (|2013[ ) 
propose a larger gas density than we assume in the current work. 
Tvib,v = 3 — 2 j Tvib,t) = 4 — 3 ? and b,u = 5 — 4 are only plotted 
up to 5R* because the lines of these bands are mostly formed 
in Region I and no information on larger distances to t h e star 
ca n be retrieved from the m. References: (1) De Beck et al.l (12 012) 
f2)lAgundez et alJ d2012h (3) ISchoier et all ( 20071) f4’) lBovle et al.l 
dl994h 


collisional rates discussed in lYelilla Pr ieto et al. (2015), 
where the calculations performed by [T obola et al.l ( 20081 ) 


for the colliding pair 28 Si 32 S:He were improved and extrapo¬ 
lated to high ro-vibrational levels. As far as we know, there is 
no data in the literature about ro-vibrational transitions due 
to collisional processes between 28 Si 32 S and H 2 . Thus, we as¬ 
sume that the collisional rates are the same for H 2 and for 
He, and that these rates, calculated for a kinetic temperature 
ranging between 300 and 1500 K, can be extrapolated up to 
~ 2300 K following the power-law T k ' . We estimate that 
the maximum 28 Si 32 S ro-vibrational collisional rates ranged 


of bands v = 1 — 0to5 — 4. Hence, n cr it ^4x 10 12 cm -3 
at the stellar photosphere, two orders of magnitude larger 
than the estimated gas density. However, the ana l ysis o f 
the HNC emission performed by [Cernicharo et all (|2013h 
is compatible with an H 2 density very close to the star 
larger than 10 12 — 10 13 cm -3 as a consequen ce of the mat¬ 
ter levitation due to t he s t ellar pulsation ( BoweiJ 
WjllacrY_fe_C^nd]^fl|_|l^^ fe Cernicharo! 


Agundez et al. 20121 ; Cherchneflll2006l l. being similar or ex- 


19881: 


2006i: 


ceeding the estimated critical density (n gas (R*) > n cr i t ). 
Thus, the vibrational excitation would be dominated by col¬ 
lisions and 28 Si 32 S might be vibrationally under or close 
to LTE in the vicinity of the stellar photosphere, mean¬ 
ing that the vibrational temperatures at this position we 
derive in the current work are significantly constrained by 


the power law consider to prevail throughout Region I and 
should be assumed as lower limits. In this case, a more re¬ 
alistic approach would include a new Region adjoining the 
photosphere of less than 1R* thick (e.g., an extension of 
the photosphere), where the dependence of the vibrational 
temperature with the distance to the star would be much 
steeper than in the rest of Region I (Fig. [4]). The larger flat¬ 
ness in Region I showed by T v ib,n = 3 - 2 , T v ib,„ = 4 - 3 , and 
T v ib,v = 5-4 in comparison to the vibrational temperature 
derived from bands v — 1 — 0 and 2 — 1 are probably re¬ 
lated to the lack of reliable information about the emission 
component of the lines of bands v = 3 — 2, 4 — 3, and 5 — 4. 
Thus, the vibrational temperature for these bands could be 
steeper th an our best fit suggests th roughout Region I, as the 
results bv lVelilla Prieto et al.1 (I2015T) indicate. On the other 
hand, the decrease experienced by the vibrational temper¬ 
ature in Region I when the vibrational quantum number v 
increases (T vi b,„ = 1 _ o{R+) compared to T v ib,u = 2 - i(R*), 
T v ib,u = 3 - 2 (R*), and so on) derived from our fits seems to 
be significant, as the uncertainties suggest (Table [2]). This 
variation is probably caused by the different A-Einstein co¬ 
efficients of the ro-vibrational transitions of bands v = 1 — 0 
to 5 — 4, that are about i) up times larger than for band 1 — 0. 


4-1-4 The red-shifted, emission excess 

The observed lines of the v = 1 — 0 band with Ji ow < 40 
show a red-shifted emission excess at Doppler velocities of 
>6 — 10 km s^ 1 with respect to the lines derived from 
the best fit. The affected lines involve ro-vibrational lev¬ 
els with energies below 700 K, suggesting that the emitting 
gas is located further than ~ 15R* from the star. This sce¬ 
nario is compatible with the lack of an emission excess in 
higher excitation lines, which are mostly produced closer to 
the star than ~ 10R*. A simi lar red-sh ifte d em ission ex¬ 
cess was already reported bv lFonfrfa et al.l (120081 ) in the ro- 
vibrational lines of C 2 H 2 and HCN. Interferometer obser¬ 
vations carried out with HPBW ~ 0.25 — 0.35 arcsec of the 
lines Si 34 S(u = 0, J = 15 — 14), H 13 CN(u = 0, J = 3-2), and 
SiO(« = 0, J = 6 — 5) (IShinnaea et al.1 [20091 : IFonfrfa et al.l 
120141) also showed red-shifted emission excesses probably 
prod uced in the dust formation zone (1 7th ^ r < 30R*; 
e.g.. IFonfrfa et akll2008l : iDecin et al .11201(1 ). the same region 
where the ro-vibrational li nes prese nted and analysed in the 
current work and bv lFonfrfa et al.l (l2008ll are formed. 

A spherically symmetric envelope model cannot repro¬ 
duce the red-shifted emission excess and the absorption com¬ 
ponent of many lines at the same time. Since these absorp¬ 
tion components highly constrain the possible structure, the 
expansion velocity held, and the physical and chemical con¬ 
ditions prevailing throughout the envelope, other models 
should be invoked to explain the observations. As most of 
the emission component of the SiS lines of the v = 1 — 0 
band comes from the rear part of the envelope close to 
the star while the absorption component is formed in front 
of the star, the red-shifted emission excess could be mod¬ 
elled assuming an asymmetric SiS abundance and/or vi¬ 
brational temperature distributions. This asymmetry could 
be a consequence of (1) gas expanding at velocities higher 
than expected in the surroundings of the stellar photosphere 
t hrou gh, e.g ., the existence of an outflow behind the star 
(Fonfrfa ct ah|201-|), (2) the influence of the companion on 









































































SiS in the innermost envelope of IRC+10216 9 


Table 4. Isotopic ratios in the innermost envelope 


Ratio 

Value 

Error 

Reference 

28 Si/ 29 Si 

17 

+ 5/-4 

1 


18 

+ 2/-2 

2 


8.2 

+ 1 . 7 / - 1.7 

3 


15.1 

+ 0 . 7 / - 0.7 

4 


17.2 

+ 1 . 1 / - 1.1 

5 


> 15.4 

— 

6 

32 S/ 34 S 

14 

+6/-4 

1 


22.0 

+2.5/ - 2.5 

2 


30 

+8/-8 

3 


19.6 

+1.3/ - 1.3 

4 


18.9 

+1.3/ - 1.3 

5 


21.8 

+2.6/ - 2.6 

6 


enough to derive the isotopic ratios 28 Si/ 29 Si and 32 S/ 34 S 
with a relatively good accuracy of ~ 30 and 40 per cent, re¬ 
spectively (Table0. The fit to the only 30 Si 32 S line suitable 
to be analysed gave an isotopic ratio 28 Si / 30 Si of 39 with an 
error of 60 per cent. The lack of available lines and the large 
uncertainty prevent us to derive reliable conclusions about 
this isotopic ratio. 

The ratios 28 Si/ 29 Si and 32 S/ 34 S are compatible with 
previous estimates derived from observations carried out in 
the mm wavelength range (Table [4j. These observations 
mostly recovered the long scale emission of the envelope 
that comes from la yers composed of matter expelled more 
than ~ 800 yr ago dCernicharo et aki 12015a ). while our ob¬ 
servations traced gas recently ejected. Hence, no significant 
changes in the nucleosynthesis of the most abundant Si and 


(1) This work (2) 

Agundez et al. 

2012 

) (3l|Patel et al.l(201lh (4) S isotopes are expected to have happened during the last 

Patel et alJ (200£ 

) (5) 

He et al. 

2008 

) (6) Kahane et all i 2000) 1600 yr. 

and Cernicharo et al.l ( 

2000|). The isotopic ratios by Patel et al. 


d201lh presented here were calculated by assuming as reliable the 


isotopic ratios derived from optically thin lines that do not show g CONCLUSIONS 

extended emission when possible (see table 6 of this reference). 


the ejected matter by the AGB star by its orbital movement 
llCernicharo et aDl201Bal lch. or (3) an anisotropic mass loss 
from the AGB star. 

A simple 2D envelope model can be pro posed and its 
mole cular emission calculated with our code (Fonfrfa et al.1 
120141) in order to show the feasibility of an asymmetric emis¬ 
sion. In this model, we defined a conical region behind the 
star with its axis parallel to the line-of-sight, its vertex 
matching up with the centre of the star, and an aperture 
of 120°, where the SiS abundance and the vibrational tem¬ 
perature for band v = 1 — 0 could be different than in the 
rest of the envelope. These magnitudes could only depend 
on the distance to the star in both regions. The aperture was 
chosen to modify the part of the emission component of the 
lines of the v = 1 — 0 band where the excess was found. The 
initial physical and chemical conditions were those derived 
from the best fit assuming the spherically symmetric enve¬ 
lope (Tables |T] and O. In the best fit adopting the asymmet¬ 
ric model (see Fig. 0), the vibrational temperature of the 
v = 1 — 0 band in the conical region behind the star was 
incremented from ~ 230 up to ~ 300 K at the outer accel¬ 
eration shell (located at 201 ?*), following the law oc r -0 ’ 45 
outwards. The variation of the SiS abundance in the conical 
region derived from the fitting process was negligible. This 
2D envelope model did not introduce any significant mod¬ 
ification in the synthetic lines of bands v = 2 — 1, 3 — 2, 
4 — 3, and 5 — 4 compared to the results of the spherically 
symmetric envelope model. 


4.2 Isotopic ratios 

The observed lines of 29 Si 32 S, 28 Si 34 S, and 30 Si 32 S belong to 
band v = 1 — 0 and show shapes compatible with P-Cygni 
profiles (Fig. |T|). These lines are weak, with differences be¬ 
tween the maximum emission and absorption of about three 
to four times a for lines with Ji ow < 50. Most of the lines in 
the observed frequency range are blended with C 2 H 2 , HCN, 
and 28 Si 32 S lines. Only a handful of lines over the detec¬ 
tion limit are sufficiently isolated and can be conveniently 
fitted (Fig. [I]). The lines of 29 Si 32 S and 28 Si 34 S were strong 


In this paper, we have presented high spectral resolution ob¬ 
servations in the mid-IR range (~ 12.5 — 14.0 /im) towards 
the AGB star IRC+10216. Two hundred and four lines com¬ 
ing from band v = 1 — 0, 2 — 1, 3 — 2, 4 — 3, and 5 — 4 of 
28 Si 32 S, 26 from v = 1 - 0 of 29 Si 32 S, 20 from v - 1 - 0 of 
28 Si 34 S, and 15 from v = 1 — 0 of 30 Si 32 S have been identi¬ 
fied in the observed spectrum. About 30 per cent of the total 
set of detected lines are unblended or partially blended and 
could be fitted and analysed. We have found that: 

• the quality of the fits is higher by assuming that the 
gas expansion velocity profile grows linearly from the photo¬ 
sphere up to 51?* following the dependence ti e xp(17?* ^ r < 
51?*) = 1 + 2.5(r/7?* — 1) km s _1 compared to the fits based 
on a constant velocity of 5 km s _1 . Our fits support the 
expansion velocity profile proposed bv lFonfrfa et alJ (120081 1 
outwards from 57?*, i.e., u e x P (57?* SC r < 207?*) = 11 km s _1 
and u e xp(r Js 207?*) = 14.5 km s _1 . 

• The abundance of 28 Si 32 S is 4.9 x 10 -6 between the stel¬ 
lar photosphere and the inner acceleration shell (located at 
57?*) and decreases linearly down to 1.6 x 10 ~ 6 at the outer 
acceleration shell (at 207?*). The observed lines are compat¬ 
ible with an abundance of 1.3 x 10 ~ 6 at 507?*. The decrease 
in the abundance beyond the inner acceleration shell is prob¬ 
ably explained by the depletion of SiS on to the dust grains 
due to its refractory nature and the lack of gas-phase chem¬ 
ical reactions expected to be enabled in this region of the 
CSE. 

• 28 Si 32 S is vibrationally out of LTE in most of the enve¬ 
lope. The vibrational temperature for the vibrational ground 
state is 1560(7?*/V) 0 ' 58 K between the stellar surface and the 
inner acceleration shell (Region I) and 613(7?*/r)° ' 2 K up 
to the outer acceleration shell (Region II). For the v = 1 
vibrational state, the relation is 1290(7?*/r) 0,60 K in Region 
I and 490(7?*/r) 0 ' 68 K in Region II. Beyond the outer accel¬ 
eration shell, the vibrational temperature for these states is 
compatible with an exponent of —1. For v = 2, 3, and 4, the 
vibrational temperature is compatible with 760(7?*/r) 0 37 , 
550(7?*/r) 0 ' 35 , and 312(7?*/r) 0 ' 20 K in Region I, respec¬ 
tively. 

• 28 Si 32 S can be assumed to be rotationally under LTE 
with a rotational temperature of 2330(7?*/r) 0 ' 62 K between 
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the stellar photosphere and the inner acceleration shell and 
860(i?*/r) 0 55 K outwards. 

• There is a red-shifted emission excess at large velocities 
in the 28 Si 32 S lines of bands v = 1 — 0 and 2 — 1 not noticed 
in the lines of bands v = 3 — 2, 4 — 3, and 5 — 4. This 
excess can be explained by an asymmetry of the vibrational 
temperature distribution for what the gas behind the star 
shows a higher excitation around ~ 20 .R* that mostly affects 
the lines of band v = 1 — 0 . 

• The isotopic ratios 28 Si/ 29 Si and 32 S/ 34 S derived from 
our observations are 17 and 14, respectively. These ratios, 
referred to matter in the vicinity of the star, are compatible 
with the previously proposed values for outer shells of the 
envelope. 
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APPENDIX A: GAS EXPANSION VELOCITY 
FIELD 

In order to fit the observed lines, we initially c hose t he ex¬ 
pansion velocity profile derived bv lFonfria et ah ([2008) from 
the analysis of several hundreds of C 2 H 2 and HCN lines in 
the mid-IR. During the current work we slightly modified 
this velocity profile to improve the quality of the fits to 
th e SiS lines provin g that the velocity profile derived by 
iFonfrfa et al.1 (120081 ) i s essentially comp atible with the ac¬ 
tual profile. However, IDecin et al. l (20141 ') has recently pub¬ 
lished a new profile based on the analysis of high angular 
resolution observations (HPBW ~ 0.33 arcsec) carried out 
with ALMA in the submillimetre wavelength range. In this 
new profile, the gas expands at 2 — 3 km s -1 between the 
stellar surface and 5.6i?* and grows linearly up to Hi?*, 
where it reaches the terminal expansion velocity (profile 1; 
not shown). These authors also proposed a more complex 
profile compatible with their observations where the expan¬ 
sion velocity is 8 km s _1 between 8 and 10i?* (profile 2; not 
shown). 

We tried to fit our observations with these new veloc¬ 
ity profiles finding that it is not possible to get better fits 
than that the obtained with the profile derived in the cur¬ 
rent work (Section [4]T7T|. The largest deviations of the syn¬ 
thetic from the observed lines appear in the lines of the 
ro-vibrational band v = 1 — 0 with Ji ow ranging from 30 to 
50 (Fig. IA1I) . The maximum absorption of these lines oc¬ 
curs between —10 and —9 km s -1 . Assuming a spherically 
symmetric expanding envelope, the maximum absorption is 
produced in front of the star by the shells of the envelope 



Figure Al. Comparison between the observed 28 Si 32 S line 1 — 
0R(38) (black histogram) and t he synthetic l ines c alculated with 
the expansion velocity profiles of lDecin et al.l (12014! ) (profile 1 and 
profile 2 in the text in solid and dashed red lines, respectively) and 
the current work (blue solid line). The synthetic lines are derived 
from the best fits to the observed lines of band v = 1 — 0. The 
curves in grey depart 3 a from the observations, where a is the rms 
of the noise of the spectrum (~ 0.15 per cent of the continuum). 


located at distances to it smaller than ~ 15 — 20J?* and 
where the overlapping fraction between the molecular ab¬ 
sorption of contiguous shells is maximum. A gradient in the 
expansion velocity field would diminish this fraction and, 
consequently, the optical depth, making the absorption less 
prominent. This diminishing could be balanced with a larger 
molecular abundance but it would also affect the shape of 
the whole absorption component. Adopting a continuous 
abundance profile described by its linearly connected val¬ 
ues at 1, 5.6, and lli?* for the profile 1 and at 1, 5.6, 8, 
10 , and lli?* for the profile 2, the best fits cannot properly 
reproduce the absorption component of the observed lines 
(Fig. lAll) . Both profiles produce absorptions for the lines of 
band v = 1 — 0 with Ji ow = 30 — 50 comprising two velocity 
contributions approximately at —13 and between —9 and 
—7 km s -1 , which are produced by the gas expanding at the 
terminal velocity and at 7 —9 km s -1 , respectively. However, 
the minima of the absorption components of the observed 
lines are between both contributions. The main characteris¬ 
tics of the absorption component of the observed lines are 
better explained with the expansion velocity field derived in 
the current work. 


























